Quantum states of a binary mixture of spinor Bose-Einstein condensates 
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We study the structure of quantum states for a binary mixture of spin-1 atomic Bose-Einstein 
condensates. In contrast to collision between identical bosons, the s-wave scattering channel between 
inter-species does not conform to a fixed symmetry. The spin-dependent Hamiltonian thus contains 
non-commuting terms, making the exact eigenstates more challenging to obtain because they now 
depend more generally on both the intra- and inter-species interactions. We discuss two limiting 
cases, where the spin-dependent Hamiltonian reduces respectively to sums of commuting operators. 
All eigenstates can then be directly constructed, and they are independent of the detailed interaction 
parameters. 
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The existence of spin degree of freedom is fundamental 
for elementary particles. Its values determine the quan- 
tum statistics of atoms as bosons or fermions. In the 
first experimental realization of a dilute weakly interact- 
ing atomic Bose-Einstein condensate (BEC), this intrin- 
sic degree of freedom is frozen by the external magnetic 
(B-) field used to spatially confine atoms. Optical traps, 
on the other hand, are capable of equal confinement of all 
spin components. Remarkable experimental progresses 
in recent years have stimulated many studies of atomic 
spinor condensates HUH- 

The simplest example of a spinor condensate is the 
spin-1 condensate [5|-|9(, whose properties are determined 
by the two symmetric spin-dependent s-wave scattering 
lengths. Two widely used atomic species are 87 Rb and 
23 Na atoms, both dominated by density dependent inter- 
actions in comparison to the spin exchange interactions. 
This prompts the single spatial mode approximation 
(SMA) for all spin components. The ground spin state 
is thus solely determined by the sign of spin exchange 
interaction parameter: being ferromagnetic ( 87 Rb) and 
polar ( 23 Na) respectively for negative and positive spin 
exchange interactions. The full quantum calculations fur- 
ther reveal paired spin singlets 0-1H, E3 ; different from 
the mean-field (MF) picture. For spin-2 condensates 
[1, [l(| 03 , three different s-wave scattering channels ex- 
ist, giving rise to two spin-dependent parameters. The 
ground state can take the cyclic phase in addition to the 
ferromagnetic and polar phases. Spin singlets remain 
possible although they can now be formed by either two 
or three atoms [1, 0, [ll], [13]. The spin ground states 
become more complicated for condensates of atoms with 
larger spins, such as the spin-3 52 Cr atoms [J]. 

A spinor condensate can display remarkable non-MF 
features [l3j . such as the anomalous quantum fluc- 
tuations of different spin components and their quantum 
entanglement. These strongly correlated states are also 
observed for quasi-spin 1/2 bosonic gases 

EMI, 

realized 

by two distinguishable atomic species, or two internal 
states of a single species conditional on the conservations 
of numbers of atoms in each spin component (l7l Il8j . 



This work concerns ongoing experimental studies of 
mixtures of spin-1 condensates. It combines the features 
of two-component bose gases and spin-1 condensates. 
Previously, the ground state phase diagrams of this type 
of mixtures have been studied by us using simulated an- 
nealing under MF approximation and SMA, comparisons 
with the quantum states from full quantum diagonaliza- 
tions found general agreement with the MF states [l9j . 

In this study we hope to understand the quantum state 
structures including spin fluctuations for a mixture of two 
spinor condensates. We adopt ^M F (r) (<&j\/ F (r)) (Mp=- 
1,0,1) as the field operator that annihilates a boson of 
species one (species two) at position r in spin component 
Mp. The interaction between two distinguishable spin-1 

( 12) 

atoms arc parameterized by the scattering lengths a y Q 1 ' 2 
in the channels of total spin F to t = 0, 1, and 2 respec- 
tively. The corresponding pseudo-potential is given by 
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[HHl, where g { % = 4rfa^"' 2 //i. /i is their reduced 
mass. 7^0,1.2 is the corresponding projection operator, 
and F\ ■ F 2 = V 2 — V\ — 2Vq. An alternative form is 



Vi 2 (ri - r 2 ) = -(a 



■jflFi-Fa+TWOn-ra), (1) 



where a = (g\ 
,(12) 



(12) 



j=(2g^>-3g[ 12) +g^)/2. 

Using the SMA with mode functions ^(r) and <f>(r) 
for the two atomic species, the field operators can be 
expanded as 



,(12) 



(-. 9l 12) + 



g { 2 12) )/2, and 



#j(r) = aiV(r), $i(r) 



(2) 



where en and hi are the respective annihilation operators 
for the spin component i. They satisfy the usual boson 
commutation relations. In the absence of B-field, the 
spin-dependent Hamiltonian for the mixture model we 
discuss becomes 
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where the interaction parameters C\, C 2 , and C\ 2 take 
the same definitions as in [19j . and the inter-species 
singlet pairing operator is defined as 8{ 2 = aJSLj — 
aJSJ + aLi&i- The atom numbers for each species are 
conserved, thus operators 7Vi i2 are constants N^' 2 \ 
We have defined operators Li_ = \/2(a\ao + aja-i), 
Li+ = >/2(a(joi + aliflo)) an d L\ z = (ajai — alifl-i); 
L 2 _ = v/2(SjS + &JS_i), L 2+ = V2(6 bx + bljo), and 
L 2z = — SLi&— 1) , which follow angular momentum 
algebra 0, [H[ . Thus Lf and L !Z have a complete set of 
common eigenvectors \U,li Z )- 

The first three terms commute with each other, while 
all of them are non-commuting with the fourth term in 
([3]). This implies the eigenstates will depend on interac- 
tion parameters, different from single spinor condensate. 
We will consider two limiting cases; First with 7 = as 
was studied before in the Ref. [22| and secondly when 
Ci/?i = C2P2 = Ci2/3/2. Both corresponds to situations 
where the spin-dependent Hamiltonian reduces to sums 
of respectively commuting operators, thus allowing for 
analytical derivations of all eigenstates. 

The case of 7 = ignores the inter-species singlet pair- 
ing interaction. The eigenstates for our model Hamilto- 
nian (|3|) arc then simply given by the simultaneous eigen- 
states of the four operators L 2 , L 2 , L 2 = (L\ + L2) 2 and 
L z = L\ z + L2 Z denoted by Even with a 

nonzero B-field, including linear Zeeman shifts, the ex- 
act eigenstates for a spin-1 condensate are known Q 

as IZi.Ii,>! = Z~ l/2 (!,!_) h ~ h * (4 2)t ) ?1 |vac), 

|M2*>2 = Z 2 - 1/2 (i 2 _)' 2 "^(6t) i2 (^ 2)t ) 92 |vac>, 

where i^ 2)t = aj 2 - 2a\al 1 , = bf - 2b[V'_ ll Z 1 

and Z2 are normalization coefficients. The total num- 
bers of atoms in each species are now constrained by 
NV' = lj + 2qj for j = 1,2. A complete set of basis 
for the mixture of two spin-1 condensates can be con- 
structed from the product state |Zi, h z )i ® |Z 2 , Z 2z ) 2 . This 
corresponds to coupling of two angular momentums and 
we find \h,h,l,l x ) = I] j la . W2 JZi, Zi 2 )i ® |Z 2 , Z 2z ) 2 

mi, 771,2 

where C,' / , , denotes the Clebsch-Gordan coefficient 

(hi hzl h, foz\l, h) 
= (-l)'-^V2T+lf Z 1 / 2 \ V (4) 

The last factor with the parentheses is 
the familiar Wigner 3j-symbol. The cor- 

responding eigenvalue is given by E = 
id/?! [h(h + 1) - 2JV«]+iC 2 /3 2 [Z 2 (Z 2 + 1) - 2iV( 2 )] + 
\C 12 f3 [1(1 + 1) - Zi(Zi + 1) - Z 2 (Z 2 + 1)], under the con- 
strains of < Zi < 7V (1) , < Z 2 < 7V (2) , and 
I Z x — Z 2 | < I < l\ + Z 2 . The ground state is easily found 
and its quantum fluctuations studied by varying the 
intra- and interspecies spin-exchange interactions as well 



as the linear Zeeman shift. These and other relevant 
details are not the main focus of this study and will be 
published elsewhere p3l |. 

For large anti-ferromagnetic spin-exchange interaction 
between the two species, the MF ground state is termed 
the AA phase corresponding to fully polarized spins of 
each species along opposite directions, irrespective of the 
natures of intra-species spin exchange interactions [l9j |. 
In the quantum treatment, the above MF ground state 

becomes Z~ 1 l 2 a}- i N |vac). It fails as it is not an 

eigensate of the Hamiltonian. With the help of 2L\ -I/ 2 = 
(L1+L9) 2 — L\ — L 2 , the ground state is found to take the 
form \N^\M 2 \ |ATM - N^\,l z ) with -|AT« - N^\ < 
lz < |ZV (1) - N (2) | when 7 = in the AA phase. For 
the special case of equal populations in two species with 
A^ 1 ' = iV^ 2 ) = N, it reduces to the strongly correlated 
form 

^° A= J2N + 1 ^ (-l) N - m \N,m)i®\N,-m) 2 ,(5) 

m=—N 

which is maximally entangled as was discovered before 
for a spin-1 condensate in a double well [24|. It appears 
in our model as the ground state, instead of a dynam- 
ically created state from adiabatically tuning the spin- 
exchange interaction from antifcrromagnetic to ferromag- 
netic through an optically induced Feshbach resonance 

Without loss of generality, we consider the more gen- 
eral case of unequal populations of > N^ 2 \ The op- 
erator L 2 is found to be lower bounded by I = N^ 1 ' — N^- 2 ' 
instead of the nominal minimum eigenvalue 1(1 + 1) = 0. 
The corresponding state changes into 

*-<-irv3+TE(££-S.) 

x|Zi,Zi 2 )i ® |Z 2 ,Z 2z ) 2 , (6) 

where h = l 2 = N (2) , Z = - N<- 2 \ and l z = 

hz + hz ( — Z < l z < I). The 3j-symbol is given by 

(it-lH- 1 '"-'"* 

(2l2)\(2l)\(h+l lz )\(h-l lz )\ 1 1/2 
_(2Zj + 1)! (Z 2 + l 2z )\ (h - hz)\ (I + h)\ (I - lz)\\ ■ 1 ' 

In the absence of a B-field, states iPaa with different l z 
are degenerate. 

The analogous state ■j/'aa f° r ^ w0 pseudo spin-1/2 Bose- 
Einstcin condensates was discussed before constructed 
from two orbitals [TtJ or two atomic species [2{|. The 
angular momentum like state then reduces to 

\ a \) \ a i) 

\h,l u ) = ^= J V JZ_^ = |vac), (8) 

V(h + hz)W(h-hzy. 
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where a-j- and are the annihilation operators for one 
condensate. The state \l 2 ,l 2z ) is constructed similarly 
but with the two boson operators being 6-j- and 64. for the 
other condensate. When these two states are substituted 
into the state of Eq. (J6j> , we arrive at a simpler form as 
that in [H El with 



ll z 

AA 



z -l/2 



l+l z 



l-l. 



a\b\ 



4 b \ 



(9) 



where we assume l\ > I2 and I = Z x — Z 2 - This simple 
yet elegant state displays strong correlation in a form 
that our state V'aa l° r ^ w0 s P m ~l condensates cannot be 
reduced to. 

We now consider the second special case with C\fi\ = 
C2P2 = Ci2/3/2 with the reduced model Hamiltonian 

h s = - A c 12 pi 2 + lc 12 jQ\ 2 e 12 - \c 12 p{Ni + N 2 ).(10) 

4 o I 

As in Ref. [1, [n|, we define operators iS_ = <Sj_ = 612 
and S z = (Ni + N 2 )/2 + 3, which satisfy the SU(1,1) 
commutation relations, 



[S X ,S±]=±S ± , [S+,S-] = -2S g . 



(11) 



The Casimir operator S 2 that commutes with S± and S z 
-S+S_+S 2 Z -S Z . (12) 



is given by @, [ll( 
S 2 



All terms in the Hamiltonian of Eq. (|10|) again are found 
to commute with each other, which makes it possible to 
find its eigenstates as the simultaneous eigenstates for 
the operators L 2 and Firstly, we consider simul- 

taneous eigenstates for operators S 2 and S z , denoted as 
|<S,«S Z ) with respective eigenvalues S(S — 1) and S z . Be- 
cause S z remains positive and S+S- = S 2 — S z — S 2 is 
positive semidefinite, the allowed values are S = no/2 + 
3 (no = 0, 1, 2, ... ) where the minimum of S is limited by 
the constrain S z > S and S z = S + n s (n s = 0, 1, 2, . . . ). 
Since min(iS z ) = 3, min(<S) = 3. The constrain from 
atom number conservation A^ 1 ) +-/V' 2 ' = 2n s +no further 
limits the values of no and n s . n s can thus be interpreted 
as the number of spin-singlet pairs between two species, 
and rig as the number of remaining bosons. 

The model Hamiltonian form Eq. (fT0|) for this second 
case is analogous to that of a spin- 2 condensate [H [111 ]. 
Thus we can find its exact eigenstates and eigenvalues in 
a similar fashion while noting the difference of one extra 
degree of freedom from the conserved numbers of atoms 
for both atomic species. As the operators S± commute 
with L 2 and L z , the eigenstates can be classified accord- 
ing the quantum numbers no, n s , Z, and l z . Thus they can 
be denoted as |no, n s , I, l z ; A), where A = 1,2, . . . ,g no ,i is 
used to label the degenerate orthonormal manifold, and 



E = ^C 12 pi{l + l) + -C 12 'yn s (n s + no + 5) 



-JCi2/3(iV (1) + JV< 2 >). 



(13) 



For a given set of {n , n s , I, l z }, the number of degenerate 
states g no ,i is independent of n s and l z as in the spin-2 
case [ll|, except now no = nio+?i2o, "10 = N^—n s , and 

"20 = 

AT(2) 

— n s . g no j can be found with its generating 
function defined as 



00 00 00 



G(x,y,z) = E E^o.W*" 10 / 120 ^ (14) 



nio— 0n2o— /— 



Following similar procedure outlined in Ref. |ll| , we ob- 
tain the generating function 



G(x,y,z) 



1 + xyz 



(1 - xz)(l - yz)(l - x 2 )(l - y 2 ) ' 



(15) 



Expanding G(x,y, z) around x — y — z = then gives 
the degeneracy factor A of the state |no, n s , I, l z \ A). 

To construct the explicit form for the eigenstates, we 
compute the generating function G g (x, y, z) for the max- 
imum spin states \l, l z = I), which is 

G g (x,y,z) 
00 00 

2_^h Nm )JV( 2) tl x y z 

ATd,2) = o 1=0 

= ^-A d( {1 7 cl) n (1 - < J ') _1 (i - yCT 1 



2-ni 



3=-l 

1 + xyz 



(1 - xz)(l - yz){l - x 2 )(l - y 2 )(l - xy) 



, (16) 



where h N{ i) N{2 ) ; is the total number of maximum spin 
states \l, l z = I) for a mixture of spin-1 condensates, and 
the loop £ is along the unit circle. By expanding the 
generating function G g (x,y, z), we find the six building 
blocks for constructing the eigenstates \l,l z = I), 



A?" 


= a[, 




3(2)t 


»+2 

= a - 


2a\a<_ 


6™ 


= b{, 








2616^ 


f(M)t 
1 


= ©L, 




f(M)t 
1 1 




'a\bl - 



(17) 

The general structure of the state \l,l) is then given by 

iz,o=E c ({M J {^},W)(4 1)t ) ,il (4 2)t )" 2 

x W (^T 1 (^'^"ivac},^) 

where Ui, V4, and iUj satisfy the following constrains 



iti + 2u 2 +wi +w 2 = N (1 \ 
v 1 + 2v 2 + w 1 +w 2 = N {2 \ 
ui + v\ + w 2 = I, 



(19) 
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and additionally w 2 = 0, 1 from the generating function 
of G g (x, y, z). For a fixed I, the states 



\ui,u 2 ,vi,v 2 ,w 1 ,w 2 ) 



1/2 Q(Dt 



-4 



(2)t 



x B, 



vac) (20) 



satisfying the Eq. (fT9j) forms a subspace 23. The con- 
struction of the eigenstates \no,n s ,l,l z ; A) follows Ref. 
[Til ]. We consider a series of subspaces Hy z =i) = D 



D 



with W spanned by states satisfying n s > j. 



A new subspace 2?' can be constructed by projecting the 
basis in the subspace B of definite values for n s . It can 
be simply realized by 

(P (Wl) - P {W1 + 1) )\ U1 ,U 2 , V U V 2 ,W U W 2 ) 
= (P (0) -P (lUl + 1) )|Ul,U2,Vl,«2 J ^l,^ 2 ) 

-(p(V - p(^)\ Ul> u 2 ,v 1 ,v 2 ,w 1 ,w 2 ) 



^(14)t 



Pf. 



__ o) ^ (i ( 2 )ty- ^ f )ty 



x B. 



6(2)t 



(i.i)t 



(21) 



where Pr na =o) = P^ — P^ 1 ' is the projection opera- 
tor onto the subspace with n s = 0, and the eigenvalue 
of the operator S+S- is precisely zero. We have thus 
constructed simultaneous eigenstates for the operators 
{S+S-, L 2 , L z } with l z = I. The eigenstates for L < I 
can be constructed by simply applying (L_ )'"'*, and are 
given by 



( X_)-(f(^)t)-p ( „ s=0) (i«t)- (l ( 2 )t)- 
x(p«t)-(^)t)-(f(^)t)- |vac)i (22) 

with u\, u 2 , Vi, v 2 , Wi = 0,1,2,..., oo, w 2 = 0,1, and 
Al = 0, 1, . . . , 21 satisfy the following relations 

n = ui + 2u 2 + vi + 2v 2 + 2w 2 , 

n s = wi, 

I = U\ + Vi + 2w 2 , 

l z = l-Al. (23) 



The corresponding eigenvalues are given by Eq. (JTi 

Before conclusion, we note that for the most general 
case where all interaction parameters are free to take any 
values, we can use the method of numerical diagonaliza- 
tion to find all eigenstates and the ground states. For pa- 
rameters near the two exactly solvable cases considered 
here, perturbation theory can be adopted to find the the 
approximate eigenstates and the ground states. With 
the generating function for the maximum spin states 
\l, l z = I) in hand, and recognizing that they form a sub- 
space B, the diagonalization can be carried out within 
each subspace B to find the eigenstates for l z = I. The 



states for l z ^ I can be generated by applying the oper- 
ator L_ to the eigenstates of l z = I. This way we gen- 
erate the complete structure of the eigenstates, despite 
of that the model Hamiltonian contains non-commuting 
operators and the actual eigenstates depend on the in- 
teraction parameters. 

In conclusion, we have discussed two special cases of 
exact quantum states for a mixture of two spin-1 atomic 
condensates: one without inter-species singlet pairing in- 
teraction and the other when C\f3\ = C 2 f3 2 = C\ 2 j3/2. 
Both cases reduce the spin-dependent Hamiltonian to 
sums of commuting operators, whose eigenstates can be 
constructed by finding the simultaneous eigenstates for 
all operators. For the first case with the inter-species 
pairing interaction absent 7 = 0, we have further com- 
pared the results to the MF approximation. In the inter- 
esting case of the AA phase, we find the exact eigenstates 
corresponds to a maximal entanglement between conden- 
sates of the two species. 
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No. 2009011002. 



Note added: This work has been ongoing for a while. 
It is urgent to complete the manuscript now because 
we find a recent submission to the archive by Yu Shi 
(arXiv:0912.2209) considered the same model system al- 
though focusing on the special case of 7 = 0. Our results 
agree especially in the form of the strongly correlated 
state of Eq. ((5|) in the AA phase for equal populations 
in the two atomic species. His major result that the 
state in Eq. (|5|) can be further expressed in a more tight 
form as V'aa = Z^^Tq 1 |vac), however, is wrong. 
To see this clearly, we consider the simple case of two 
atoms in each species. The AA phase now takes the 



form< A = -J 7 = r 



,(l,l)t2 



1 i(2)f fi(2)f 
3^0 n 



not of the form 



2V5 r 

given in arXiv:0912.2209. For the more general case, with 
equal populations in the two species, we can find all pos- 
sible solutions of Eq. (H3) with / = \N^ - N^\ = 
and l z = I = which forms a subspace B constructed by 
the allowed basis states \ui,u 2 ,v\,v 2 ,w\,w 2 ) . We can 
compute the matrix elements of the operators L\ and L 2 
within the subspace B, and diagonalize to find the eigen- 
states and eigenvalues with {L\) = N^(N^ + 1), and 
(L? 2 ) = (N^ + 1). The state of the AA phase in Eq. 

when 7 = is then constructed by a linear superpo- 
sition of all possible basis states in the subspace B. With 
the increasing of atom numbers, more basis states from 
the Eq. (fT9"]l will be included, and the AA phase in Eq. 
([5]) will deviate more from the state Z- 1 /2f( 1 ' 1)tw |vac). 
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